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A dynamical lattice model is used to study the viscosity and the velocity-velocity autocorrelation 
function in a microemulsion phase. We find evidence of anomalous viscosities in these phases 
(relative to water-rich and/or oil-rich phases), in qualitative agreement with other results. We 
also investigate the dynamic relaxation in the microemulsion phase. It has been suggested that 
the temporal relaxation in the microemulsion phase may be described by a stretched exponential 
Kolrausch- Williams- Watts law (</>(t) — ylexp[— (i/r)'']). In our model, we find the velocity-velocity 
autocorrelation function {v{Q)v(t)) fits this law, showing both enhanced {b > 1) and inhibited {b < 1) 
diffusion. 

PACS numbers; 68.10.Et, 66.20.+d, 64.60.Ht, 05.50.+q 



I. INTRODUCTION 

Using models to study the microemulsion phase has 
been made difScult by the combination of short-ranged 
forces, short-range order, and long-ranged disorder ex- 
hibited. A large system is needed in order to observe the 
unique isotropic "structure" , while computer power dic- 
tates what is actually feasible. Very large systems can be 
handled by Monte Carlo (MC) simulations (either con- 
tinuous or discretized), but the use of time-independent, 
non-physical moves prevents calculation of temporal ob- 
servables. Molecular dynamics (MD) allows the measure- 
ment of interesting time-dependent properties, but also 
requires large numbers of independent particles. The up- 
per capacity of MD is being continually stretched [^]-|j , 
but very large bulk systems are rarely effectively mod- 
eled. 

There have been a number approaches have been used 
to balance these concerns [pj-p^. In this paper, a dy- 
namic lattice model introduced earlier p3| is further ex- 
plored in which a Boltzmann velocity field is mapped 
onto a nearest-neighbor lattice model |^,|5[, and the 
velocity field then is integrated into the metropolis MC 
scheme. It enables the determination of properties usu- 
ally only accessible by Molecular Dynamics methods, and 
has been used to determine the diffusion constants and 
the dynamic surface tension between competing phases, 
with evidence of a corresponding wetting transition mm . 



II. DYNAMIC MONTE CARLO 

One of the aims in investigating complex surfactant 
systems is to generate a method which can be used to 
study both equilibrium and dynamic properties. Towards 
this goal, a general dynamic lattice MC model was con- 
structed by combining Monte Carlo and Molecular Dy- 
namic methods. This so called Dynamic Monte Carlo 



nonequilibrium studies. This is used to investigate the 
dynamics of a lattice Monte Carlo, while keeping the 
equilibrium states of the model intact. 



A. Description of the DMC Model 

The DMC simulation algorithm incorporates pseudo 
velocities into Monte Carlo calculations within a three- 
dimensional lattice model. In the base lattice model, cu- 
bic site i is occupied by exactly one type of the three 
chemical species in a complex surfactant system. These 
are denoted by a statistical variable cji, which can take 
on the values { — 1,0, +1}, corresponding to cubes of oil, 
surfactant, or water respectively. A vector variable Si, 
defined at each lattice site, specifies the orientation of 
amphiphile particles. The water and oil sites give no ori- 
entation. This may be accomplished by the construction 
of a hybrid variable: 
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Thus, only surfactant sites have a non-zero s.^. The vari- 
able Si is normalized as unit vector quantity that varies 
continuously in three dimensions. 

The principle advantage of lattice models, compared 
to continuum models, is that they dramatically simplify 
the simulated motions of the particles. The particles may 
be located only at lattice sites, and there are no vacan- 
cies in the lattice. Therefore, each site cannot exchange 
with other sites independently. Only pairs of sites, to- 
gether, can be considered for exchanges. This exchange 
is achieved by considering both the relative velocity be- 
tween two nearest-neighbor sites, and the potential en- 
ergy through the Metropolis algorithm. Although the 
simulation does not reflect realistic molecular sizes or in- 
dividual molecular movement, the simplicity and short 
interaction range make the lattice model far easier to 
study in detail than continuum models, and the qualita- 
tive results obtained from the model well represent the 
macroscopic behavior of oil-water-surfactant systems. 



In order to depict dynamics with this equihbrium 
model, a velocity field is apphed to the lattice MC. The 
temperature of the simulation box fixes the kinetic energy 
of the particles in the system (a form of equipartition): 
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At some time t = (really the first DMC step), the sys- 
tem is heated by randomly assigning velocities chosen 
from a Maxwell-Boltzmann distribution. The Maxwell 
distribution in one dimension is a special case of the 
Gauss(ian) normal distribution. Therefore, in each of 
the coordinate directions, the velocity components v^, 
Vy, and Vz obey the Gaussian distribution law. These 
velocities must be identified as an average or center-of- 
mass velocity of all the molecules within the lattice site, 
because of the multi-particle nature of each site. Thus, 
for a given temperature, we have 
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where a is one of the three cartesian directions. 

In the simulation, the movement of particles is achieved 
through the exchange of a pair of nearest-neighbors. As 
stated earlier, the energy of the system is independent 
of the orientation of the water and oil sites. Thus, in- 
cluding all desired features (treating water and oil sites 
as inversely symmetric, and treating the two ends of sur- 
factant molecules symmetrically), the nearest- neighbor 
Hamiltonian may be written as 
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where An = n — n' denotes the lattice difference vector 
between two sites n and n'; a and cu are the effective 
chemical potentials of water and oil; and the coupling 
constants Ca determine the bonding energies (ci = 5/2, 
C2 = 3/4, C3 = —3/8, and C4 = —1/6). This creates 
a phase-diagram which includes a microemulsion phase. 
This nearest-neighbor sum runs over pairs of sites; each 
site has 2d neighbors in d dimensions. By summing over 
all the particles in the system, the potential energy of the 
entire system is obtained. 

At each simulation step, a pseudo-random walk is cho- 
sen through simulation space. A switch variable is used 
to choose one direction {x, y, or z), which ensures that 
each pair of nearest-neighbors are visited evenly. For 
each trial exchange, the change of total energy AE is ob- 
tained by evaluating both potential and kinetic energy 
differences. The potential energy factor is the difference 
of potential energy before and after the trial exchange. 



The kinetic energy is measured by the relative velocity 
V = Vn — Wn' only in the exchange direction. This means 
only the component of the velocity in the exchange di- 
rection is considered, and the components in the other 
two directions are unchanged through the trial exchange 
(and therefore any accepted exchange). The actual rules 
give results for the two possible signs of v. for the case 
w < 0, no exchange occurs even for a favorable AE, as 
a negative relative velocity corresponds to two particles 
moving away from each other. For w > 0, a resulting 
(new) relative velocity v' is obtained: 



exp[—mv /2kBT] — 1 — exp[- 



■/2kBT] (5) 



In this manner, the final velocity components in the ex- 
change direction are determined. 

Both relative velocities and change of potential energy 
are considered for the change of total energy, and the 
acceptance or rejection of the exchange is then deter- 
mined by the Metropolis algorithm. This implies that 
the kinetic energy can supply potential energy required 
exchange, or there can be mixing between the kinetic 
and potential energy. Since the velocities are initially as- 
signed from a Maxwell-Boltzmann distribution, and re- 
tained throughout the trajectory, the potential energy 
should also be contained within a Boltzmann distribu- 
tion of all the configurations obtained. For an accepted 
trial move, the chemical species within the two sites are 
exchanged, and the velocity components in the exchange 
direction will have new values resulting from Eq. (M); for 
a rejected trial move, an elastic collision is applied, so 
that the two particles move away from each other with 
an exchanged velocity component in the direction of col- 
lision. 

The continuous spin variable Sn needs special treat- 
ment. Because surfactant molecules may move to neigh- 
bor sites, the orientation of newly created surfactant site 
must come from an equilibrated distribution to main- 
tain detailed balance in the system. This is accomplished 
by creating a "ghost" vector orientation at each lattice 
site during the relaxation procedure. The ghost orienta- 
tion at each site is also produced by Eq. (||), and accep- 
tance comes by including the potential energy through 
the Metropolis algorithm. Therefore, the ghost orienta- 
tions maintain equilibrium. When a surfactant moves 
to its neighbor site, it takes the ghost orientation of the 
site. System size is maintained by periodic boundary 
conditions. 



B. DMC Simulation in Practice 

The microemulsion phase is studied by Dynamical 
Monte Carlo (DMC) simulation. Most simulations were 
performed on systems of 24 x 24 x 24 lattice sites. A test 
of validity of the DMC method is examining the velocity 
distribution over time (DMC steps). The initial veloc- 
ities of particles are generated according to a Maxwell- 



Boltzmann distribution. At the end of the simulation, 
velocities maintain the Maxwell-Boltzmann distribution 
as predicted. Additionally, there is the requirement that 
the DMC model maintains the equilibrium states from 
the standard MC. Results for the energies of systems of 
identical compositions from the two models, lattice MC 
and DMC, were indistinguishable. 

Experiments indicate that the scattering intensity of 
microemulsion exhibits a peak at non-zero wavevector 
|J17|"P1| . This scattering behavior is considered the signa- 
ture of microemulsion, in contrast to normal disordered 
phases (which show an exponential decay with increasing 
wavevector). In the DMC, the scattering intensity is cal- 
culated as a discrete Fourier transformation of the water- 
water correlation function. We have chosen as a "typi- 
cal" microemulsion system: 26% surfactant sites, 37% oil 
sites, and 37% water sites. An equilibrated microemul- 
sion phase is obtained after 4000 to 6000 simulation full 
lattice passes (DMC steps). Simulated scattering data 
from such a phase is presented in Fig. nl 
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FIG. 1. Water-water structure function of microemulsion 
phase (T=1.5). 



III. RESULTS AND DISCUSSION 

A. Viscosity of Microemulsion 

The investigation of the viscosity of surfactant systems 
is of great interest for both scientific and industrial ap- 
plications. Surfactant solutions show complicated rheo- 
logical behavior. It has been shown experimentally that 
the viscosity of such systems is affected by many fac- 
tors such as concentration, molecular volume, presence 
of electric charges, hydrodynamic interactions, and tem- 
perature. In practical terms, the viscosity of surfactant 
solutions is being scrutinized by industry. For example, 
the flow rate is strongly affected when these solutions 
are pumped through pipes or stirred in large contain- 
ers. The viscosity is also relevant to the design of surfac- 
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tant solution-based flooding processes, which can be uti- 
lized in oil recovery (following the primary and secondary 
stages) . This is of great interest since some two-thirds of 
the oil originally in place can remain after the secondary 
recovery process. 

The DMC has been used to investigate viscosity in mi- 
croemulsion. A (unit) velocity gradient was applied to 
the system in one direction of the 3D lattice, in order to 
simulate the behavior of fluids under shear stress. The 
velocity component in x-direction (v^) is forced to in- 
crease with increasing y. Across each lattice spacing, Vx 
is increased by a constant amount (the shear rate) in the 
y-direction. The shear rate is a relative unit based on 



{vl) = l {a^x,y,z) 



(6) 



Since all the lattice sites on the same xz-j>\a,ne suf- 
fer the same magnitude of shear, the exchanges between 
nearest-neighbors in the xz-plane are allowed to be per- 
formed normally. The exchange of particles across the 
2/-direction, which exchanges between two adjacent xz- 
planes, needs to consider the different shear rates of the 
two lattice sites. The exchange of two particles by ran- 
dom motion may happen between sites in planes yi and 
2/2 (see Fig. ||). Since the particle in j/i has a lesser veloc- 




FIG. 2. Schematic of the velocities of adjacent planes in 
the DMC. 

ity component in the a;-direction than does the particle 
in 1/2, this exchange process involves a net transfer of mo- 
mentum in the direction of y. For velocity components 
in the y- and z-direction (vy and Vz), the values are sim- 
ply exchanged without any added effect. For the velocity 
component in the flow direction {vx), the relative velocity 
is calculated as before, but the net transfer of momen- 
tum between the two planes is taken into account. The 
particle that moves from plane yi to y2 is coming from 
a "slower" velocity field compared to particles in plane 
2/2- Therefore, its local velocity should be lower than its 
neighbors. On the other hand, the particle from plane y2 
should have a higher velocity than other particles in plane 
2/1 . This can be summarized by the following expression 
for the final velocity component in the x-direction: 
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where u is the shear rate. 

Viscosity has the generalized form pS 
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which uses the velocity-velocity correlation function 
Cvv(t). This can be re-written explicitly for Cyy(t) with 
the result 
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Use of the correlation function for viscosity measure- 
ments requires long lists of velocities for each particle at 
each DMC step. This large physical memory requirement 
is quite restrictive. Another way is to take advantage of 
the relative velocities used in the exchanges. Labeling 
this relative exchange velocity as w^.q, we can use the 
sum 
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This sum is over only accepted exchanges of two nearest- 
neighbors in the y-direction. L is the linear lattice size 
(generally, L = 24: in this work). The variable r is in- 
troduced as the reduced time step, which ensures that 
the average velocity is equal to the average rate of ex- 
change of the sites ((|v|)t = {P{v))). It is a constant 
of the simulation method. One of the advantages of this 
method is that it can work with periodic boundary con- 
ditions, since the shear rate of one site is only relative 
to its neighbors, and can only be considered during the 
exchange. Furthermore, since the flow is not explicitly 
introduced by increasing the velocity of each lattice site, 
and the kinetic energy is defined through local velocities 
(which are balanced during each exchange), any heating 
affect of the process is minimized. 

Results for these viscosity measurements are presented 
in Table |. These results are for an equilibrated mi- 
croemulsion phase (as determined by S{k) data), at a 
constant temperature. The temperature can be indepen- 
dently determined through the kinetic energy. One can 
see that the apparent temperature falls as the shear rate 
is raised. This is due to the increased transference of 
kinetic energy to potential energy, also reflected in a cor- 
responding increase in the potential energy. 

TABLE I. The viscosity from DMC of the microemulsion 
phase at T = 1.0. 



Shear Rate 


Shear Viscosity (r;) 


Internal Energy 


Final r 




((mJ)i''V^) 


{J) 


(J/fcB) 


0.005 


2.282 


-1.448 


0.979 


0.010 


1.862 


-1.450 


0.983 


0.015 


1.523 


-1.461 


0.983 


0.020 


1.371 


-1.465 


0.974 


0.025 


1.305 


-1.469 


0.969 


0.030 


1.227 


-1.474 


0.965 



DMC calculations of the shear viscosity of microemul- 
sion with shear rate changes exhibit some notable fea- 
tures, illustrated in Fig. 0. This bears some discussion 



about what is expected for such fluids. Currently, sev- 
eral studies of viscosity of microemulsion have been pub- 
lished in the literature. Gogarty and Dreher |2^] have 
reported experimental results of microemulsion which 
shows strong non-Newtonian behavior. Many researchers 
have given values of microemulsion viscosity at only a 
single value of shear rate, or shear stress p4|-|28[|. In 
some cases, non-Newtonian or anomalous behaviors were 
observed. 124-26 28] However, there are several other re- 
searches which report Newtonian viscosities p9|-|3l|| . This 
may result from the contradictory elements of small in- 
teraction range and weak velocity variation. The hydro- 
dynamic velocity gradients applied in these experiments 
are not high enough to create a strong velocity variation. 
Thus the system needs a relatively longer time to allow 
the velocity variation to be transferred compared to the 
renewal time of the system structure. This limitation is 
not always apparent in the microemulsion phase contain- 
ing comparable amounts of water and oil. In particular, 
in a viscosity measurement of microemulsion obtained by 
a cone-and-plate rheometer, both Newtonian and non- 
Newtonian behaviors have been observed |p3]. 

The data obtained from the DMC simulation exhibit 
a strong non-Newtonian behavior, which corresponds to 
only some of the experimental results. There is still much 
work to be done in this area before one can draw any con- 
clusions as to which experimental results are correct (per- 
haps both sets are). However, we do know this behavior 
is not an artifact of the DMC, as water-rich (oil-rich) 
phases exhibit only Newtonian viscosities. On the other 
hand, there are some shear-induced perturbations in the 
simulation which may affect the viscosity results. The 
change of velocity component resulting from the shear 
is not an actual Boltzmann factor. This should affect 
the results of Metropolis algorithm in determining the 
acceptance or rejection of the exchange in the following 
simulation steps. In other words, a change of configura- 
tion is induced in the system corresponding to the shear 
rate. Evidence of this is the change of the internal en- 
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FIG. 3. The dependence of viscosity measurements on 
shear rate. 



ergy and temperature of the system which are observed 
at high shear rate. Compared to the change of viscos- 
ity with shear rate, however, the affect of shear-induced 
error does not seem considerable. At lower shear rates, 
the results are seemly linear. In any case, in order to 
minimize the perturbation of shear on the system, and 
statistical error in calculation, the viscosity is determined 
at a relatively low shear rate. Thus, we shall now fix the 
shear rate at 0.01 (relative velocity units). 

Purely thermal effects are also interesting to consider. 
Heuristically, viscosity should decrease with increasing 
temperature, but the overall behavior of microemulsion 
viscosity with temperature is not well understood. To 
investigate phase behavior of rj, the shear viscosity of 
microemulsion phase was calculated and compared with 
the oil-rich phase at different temperatures. These results 
are presented in Table || and Figure Q. 

TABLE II. Shear viscosity of the microemulsion phase and 
the oil-rich phase (u — 0.01). 



Temperature 


r] (microemulsion) 


V (oil) 






{(mjf/e) 


((mJ)i/V^) 


% Change 


1.0 


1.862 


1.979 


5.9 


1.2 


1.764 


1.920 


8.1 


1.5 


1.761 


1.850 


4.8 


1.7 


1.738 


1.812 


4.1 


2.0 


1.689 


1.764 


4.2 


2.2 


1.652 


1.736 


4.8 


2.5 


1.650 


1.699 


2.9 


2.7 


1.627 


1.678 


3.0 


3.0 
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responds to the general experimental fact that the mi- 
croemulsion phase exhibits low viscosity compared to 
other bulk phases. 

It has been found experimentally that the viscosity of 
microemulsion is strongly affected by its microstructure. 
Papaioannou and coworkers |^^ have published viscosity 
measurements of microemulsion as a function of salinity. 
They found that the viscosity of microemulsion exhibits 
two peaks: one at the oil-rich, and one at the water-rich 
phase boundaries. Additionally, a minimum was mea- 
sured in the Winsor III state. It was also found through 
other experiments that the viscosity of microemulsion de- 
creases as the structure becomes more and more bicon- 
tinuous. pj] In the DMC simulation, the microemulsion 
phase has comparable amounts of oil and water, corre- 
sponding to the bicontinuous Winsor III state. The DMC 
therefore shows good agreement with the experimental 
results. 



B. Long-time Behavior of the Velocity 
Auto-correlation Function 

Use of the time-correlation function is one of the im- 
portant approaches for describing transport processes 
and other time-dependent phenomena. This general 
function is considered to play a similar role in non- 
equilibrium Statistical Mechanics as the partition func- 
tion plays in equilibrium situations. One advantage of 
the time-correlation function is that the resulting trans- 
port coefficients are quite general and independent of the 
details of any particular model. In particular, the self- 
diffusion coefficient can be expressed in terms of the ve- 
locity auto-correlation function as 



D 



(u(0) ■ v{t)) dt 



(11) 



FIG. 4. The dependence of viscosity measurements on tem- 
perature. 



This expression is valid for any classical system in 
which the diffusion is governed by the diffusion equa- 
tion (Pick's law). Therefore, the velocity auto-correlation 
function is a valuable approach to investigate the diffu- 
sion behavior of complex systems. Autocorrelation ex- 
pressions like Eq. (^ may be nicely related to experi- 
ment. Observation of relaxation is a common experimen- 
tal method to investigate the structure and dynamics of 
systems. Because both equilibrium and dynamical prop- 
erties of the symmetric microemulsion phase are still very 
uncertain, the relaxation are of continuing importance. 

The general form of velocity auto-correlation function 
is 



One can see that with an increase in temperature, the 
viscosity of both phases decrease, which is characteristic 
of all liquids. At a given temperature, the microemul- 
sion phase shows a considerably lower viscosity than the 
oil-rich phase, and this difference decreases with increas- 
ing temperature. Such behavior of microemulsion cor- 



C™(t) = (w(0) •«(*)) 



(12) 



When t — 0, Cyy{t) simply has the value of equilibrium 
average of v'^{0) {^k^T/m by equipartition) . As time 
advances, v{t) will be less correlated with its initial value. 
Thus, the velocity auto-correlation function is expected 



to decay to zero from its initial value. In general, it is 
assumed that any C(i) decays exponentially with a time 
constant t. This can be expressed as 



C{t) = Aexp{-i/T} 



(13) 



However, it has been found experimentally that time- 
relaxation in many disordered systems does not follow 
the simple exponential decay. The microemulsion phase 
is one example of this. Losada and Lopez-Quintela Pq ] 
found a relaxation mode in microemulsion which fol- 
lows a Kohlrausch- Williams- Watts (KWW) stretched- 
exponential law. The general form of this law is 

$(i)=Aexp{-(i/r)''} , (0<6<1) (14) 

The super-exponential factor b was found to be depen- 
dent on many conditions, such as the spatial, temporal 
and energetic disorder of the system. There has been 
only limited discussion of the applicability of this law, 
and its ramifications ||3q-|3q]. 

In order to measure the velocity auto-correlation func- 
tion in the DMC model, the derivative of the mean- 
squared displacement is used: 



d{[r{t)-rm') 



dt 



{v{t)-vm 



(15) 



In the DMC method, since the movement of particles 
is achieved only by the exchange of a pair of nearest- 
neighbors and is not explicitly decided by velocity of 
each single lattice site, a "measured" {v{t) • w(0)) will 
not reflect the actual diffusive behavior of the system. 
The velocity auto-correlation function is therefore de- 
fined through Eq. (p^). 

Because of the average assembling nature of the veloc- 
ity auto-correlation function, the long-time behavior is 
more difficult to observe and less accurate Q. In de- 
termining whether relaxation follows a KWW-type law, 
the longest-time correlations are the most important. 
This requires extremely long simulation (cpu) times and 
proper equilibration before results may be examined. For 
a typical simulation, he velocity auto-correlation function 
is obtained from 3600 to 4000 equilibrated simulation lat- 
tice passes. Data obtained at different temperatures were 
fit to Eq. (Ill) for the longest part (last 2/3) of simula- 
tion time — 1000 to 4000 DMC steps (Fig. |). The fit 
is satisfactory for all microemulsion data obtained, and 
better than a simple exponential fit (Eq. (p^). Although 
all data (from t = 0) were fit satisfactorially with cither 
Eq. (|l^) or Eq. ( p^ (with 6 « 1), when only the longest- 
time data were included, the simple exponential fit was 
inadequate. Thus, Eq. (|lj) fits the entire data range 
well, while Eq. (|l^) breaks down at longer "times" . Ad- 
ditionally, auto-correlation data obtained from the DMC 
in the water-rich or oil-rich phases can not be found to 
fit a "pure" {b ^ 1) KWW law. The values of the KWW 



fit parameters are presented in Table HI. 

Both the stretched exponent b and the macroscopic 
relaxation time r are found to increase with increasing 
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FIG. 5. KWW stretched-exponential law fit to data 
of velocity auto-correlation function of microemulsion phase 
(T=1.5). 

TABLE III. Temperature dependence of the KWW equa- 
tion parameters in the microemulsion phase (A=0.1). 





A = 0.102 


T 


X = 67.2 




b = 0.577 








Temperature 


T 


h 


2.5 


145.63 


1.365 


2.2 


121.67 


1.054 


2.0 


105.79 


0.840 


1.7 


82.48 


0.630 


1.5 


69.67 


0.575 


1.3 


58.55 


0.490 


1.2 


53.65 


0.470 


1.1 


50.48 


0.459 


1.0 


48.71 


0.418 



temperature within the microemulsion phase. The in- 
crease in relaxation time may be understood as a slowing 
of the diffusion rate as the system approaches the crit- 
ical temperature. A number of reports have been pub- 
lished on the experimental study of critical phenomena 
of surfactant systems |3J]. However, the critical behavior 
in surfactant solutions is still far from being understood 
completely. One common feature found in most of the 
investigations is that phase separation occurs as temper- 
ature is raised. This may be related to the slowing of 
diffusive behavior, observed in the present model. 

The wide variation of b is of great interest. Notably, 
a value of b greater than 1 was obtained at high tem- 
peratures, which has not been observed in normal disor- 
dered phases. This corresponds to "enhanced diffusion" 
— faster than expected. Losada M] has reported a relax- 
ation study of microemulsion using a pressure-jump tech- 
nique. In this experiment, the relaxation behavior with 
6 > 1 was observed in a microemulsion phase that had 
similar concentrations of oil and water. This has been 
the only experimentally-observed example of enhanced 
diffusion in microemulsion. It is hoped that additional 
work in this area will be done to confirm this interesting 
and important result. 

The stretched exponent b is strongly temperature- 
dependent. It has been treated by the theory of criti- 



cal phenomena, and Losada found that it fit the critical 
equation 



b = A[{T,-T)/T,]- 



(16) 



The temperature dependence of b was measured by DMC 
simulation in order to confirm this behavior. The DMC 
data was fit to Eq. |6| (Fig. ^. This fit resulted in 
the value A = 0.30, and gave a critical temperature 
Tc = 3.20. From Fig g one can see that 6=1 when 
T = 2.21; at this point, the relaxation (as well as the dif- 
fusion) is "normal" and the KWW equation reduces to 
the exponential form. At temperatures below this point, 
the microemulsion shows an inhibited diffusive behavior. 
This can be related to the spatial and temporal disorder 
present in the system. It has been found in both exper- 
imental and theoretical research that diffusion behavior 
moves from anomalous to normal as the temperature is 
increased [40 4^. This corresponds to DMC results ob- 
tained for T < 2.21. As the temperature is increased 
past this point, there is an enhancement of diffusion in 
the system. 
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IV. CONCLUSIONS 

The complex systems formed between oil, water 
and surfactant molecules were examined with a three- 
dimensional lattice model enhanced by the addition of 
a coarse-grained velocity field. The shear viscosity of 
microemulsion phase was measured as a function of 
shear rate, and showed a decrease as the shear rate in- 
creased. The measurements were also performed on oil- 
rich phases, which showed typical behavior of a Newto- 
nian (classical) fluid. Compared to the oil-rich phase, the 
viscosity of microemulsion shows strong non-Newtonian 
behavior. This corresponds to at least one experimental 
observation of the existence of non-Newtonian shear be- 
havior in microemulsion [B3,E8| The shear viscosities were 



also measured as a function of temperature. In all the 
systems simulated here, the microemulsion phase shows 
a notably lower viscosity than the oil-rich phase. As the 
temperature approaches to the critical temperature, the 
difference decreases. We would expect the two viscosities 
to equal at the thermodynamic limit. These viscosity 
measurements in the microemulsion pha se qualitatively 
reproduce experimental observations |33|j34[| . 

Measurements of the velocity autocorrelation function 
were obtained by the distance displacement method. The 
most interesting observation is that the long-time tails 
of microemulsion do not obey a simple exponential de- 
cay, but rather a Kohlrausch- Williams- Watts stretched- 
exponential law. The same behavior has been observed 
experimentally in a few other disordered phases P5| , [42| . 
Microemulsion, a particular disordered phase, shows 
many properties distinct from normal disordered phases, 
which makes the long-time behavior very compelling. 
Both the stretched exponent and the macroscopic re- 
laxation time were found to be strongly temperature- 
dependent. The temperature dependence of the super- 
exponential 6, shows critical scaling. The values of h can 
be divided into three temporal ranges h = 1,6< 1,6> 1. 
Since the velocity autocorrelation function is directly re- 
lated to the diffusion constant, this suggests the exis- 
tence of normal, inhibited and enhanced diffusion in this 
model. These results provide evidence that DMC method 
employed in this work should be capable of simulate not 
only the equilibrium but the dynamics of the system as 
well. The diffusion constants and the dynamic surface 
tension between microemulsion and oil-rich phases has 
already been calculated [Q, and the measurement of 
other properties should also be possible through use of 
the DMC model. 

Due to the fact that the molecules are defined at the 
lattice sites, the movements of particles are restricted to 
unit length. Consequently, the velocities involved in the 
calculations are not completely realistic. In particular, a 
small spatial region of the velocity autocorrelation func- 
tion is anti-correlated and inaccurate at certain degree. 
Furthermore, we make no attempt to make quantitative 
comparisons of shear viscosity with experimental data. 
Such comparison would require more realistic modeling 
of the molecular interactions. 

One aspect of future work may connect velocity auto- 
correlation function results to the cooperative relaxation 
behavior of the system. The dipole correlation function 
is one way to describe of the mechanism of cooperative 
relaxation through experiment, so that the velocity au- 
tocorrelation function can be related to the relaxation 
of the dipole moment in droplets, as well as the struc- 
tural and kinetic properties of the system. The observed 
non-Newtonian behavior is a macroscopic representation 
of the microstructure of microemulsion. However, at 
present there is little theory to connect these viscos- 
ity observations with microstructural evolution. Since 
both Newtonian and non-Newtonian regimes have been 
observed experimentally, this phase progression of mi- 



croemulsion is of interest to explore further. 
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